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Abstract 

Intense experimental and theoretical efforts have recently been brought to bear on 
the determination of the symmetry of superconducting order parameter in high T c 
Cuprates as this would shed light on the nature of microscopic interactions that are 
responsible for superconductivity in these unusual systems. We use a two band model 
that reproduces the Fermi surface seen in ARPES experiments of Shen et al.,Jl| and 
study the symmetry of order parameter in the YBCO system. The model incorporates 
the effects of tunnelling in the c-direction between the planes and chains of these 
Cuprates. It is observed that a suitable choice of the phase of the pair wave- function 
in the planes and chains lead to both s- and d-wave like features. We calculate the 
detailed phase diagram of this system, observe the shift in the position of the node(s) 
of the order parameter (s) on the Fermi surface with doping and explain the seemingly 
contradictory experimental observations about the presence of both s- and/or d-wave 
symmetry [Q] in these systems. 

Quite a large number of experiments backed by careful calculations have been brought 
to bear upon the resolution of the seemingly irreconcilable and contradictory observations 
about the nature of symmetry of the superconducting order parameter (OP) in the Cuprate 
superconductors 0. An understanding of this symmetry, which is still unresolved, will shed 
considerable light on the nature of the pairing mechanism(s) in these unusual systems. 

Broadly, there are several classes of experiments that probe the symmetry of OP in the 
cuprates: i) Transport and thermodynamic measurements 0, notably the temeperature de- 
pendence of the London penetration depth (A) |4j] that show a linear temperature dependence 
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of A, indicating that at low energies the density of superconducting states (SDOS) varies lin- 
early with energy in pure samples, a signature of gapless excitations, ii) Angle resolved 
photoemission experiments that show a node (with a small offset) along the (n, it) direc- 
tion of the square Fermi surface (FS)of YE^CusCV iii) Josephson measurements: magnetic 
field dependence of dc SQUID || [/J predict a n phase shift consistent with d-wave scenario. 
Conversely, c-axis Josephson tunneling f| between twinned and untwinned YBa2Cus07 and 
Pb junction finds a Josephson current, albeit small, indicative of a symmetric OP. Recent 
atomic forced microscopy measurements of Kitazawa and his group [Q found evidence of gap 
everywhere on the FS in Bi-2212 materials. Critical current between a hexagonal grain and 
its surroundings is found to be non-zero ||10|| and it scales with the number of sides engaged, 
consistent with the presence of a conventional s-wave component, iv) Sensitivity of super- 
conducting T c with elastic scattering is found to be too weak (follows, instead, the Ioffe-Regel 
criterion, kfl — 1, where I is the mean-free path of the electrons) to conform to the behaviour 
of a pure d-wave symmetry [ PH . 

While there is a general concensus that there are nodes on the FS in all the hole- 
doped Cuprates, there is, however, no consensus about its origin. Models based on spin 
fluctuation [12| predict a d x i_ y i symmetry for the OP, while an interlayer tunneling-enhanced 



conventional superconductivity scenario |]TB| envisages extremely anisotropic s-wave gap func- 
tion. It is also useful to note that the symmetry of YBCO system is orthorhombic and not 
tetragonal, so an admixture of s- and d-components is inescapable [13], and the nomenclature 
is based on the dominance of one component over the other. 

In view of this conflicting scenario it is useful to look for models that incorporate i) 
the essence of the band structure of the system concerned and ii) reproduce much of the 
experimentally observed features. Comberscot and Leyronasfl5| have tried to incorporate 
the fact that in YBCO systems, the existence of both chains and planes can have a strong 
bearing on the nature of pairing wave function that occurs. We use a similar model and 
show that it is indeed possible to reproduce many of the experimentally observed features 
of the OP without explicitly considering a d-wave pairing function and features consistent 
with both d- and s-wave are obtainable with a suitable choice of the relative phase of two 
isotropic pairing functions residing mainly on the plane and chain derived bands. 

We consider the Hamiltonian: H = Hq + Hj where, 

H = Z^lcCLck* + t'kyAyjkya H" ^^0^ + fl.C.)] (l) 

k<7 
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6k and e' k are energy dispersions, and c^, d ka are destruction operators for electrons in the 
plane and the chain respectively. The off-diagonal term t is small (t ~ 50 meV) PE| , momen- 
tum independent and introduces band interaction. This implies that in YBCO system, an 
electron is never in a true eigenstate of planes or chains, rather, it has components on both 
the bands. This term shows up in the band-crossing region leading to anticrossing features 
with two sheets of the FS as seen in YBCO. Thus moving along a given sheet of the FS, 
an electron changes its character from being in the plane band to being in the chain band. 
This feature has its bearing on the relative sign of the OP as we discuss below. 

In order to fit our band structure with a realistic one obtained from the ARPES data|| 
for the YBCO system, we take five nearest neighbour hopping parameters in the plane, in 
both a and h direction and two nearest neighbours in the chain. The dispersion relation in 
the plane is now 



6k = — fi p — 2tiCos{k x ) + Cos(ky) + 4:t2Cos(k x )Cos(k y ) 
- 2t 3 Cos(2k x ) + Cos(2k y ) - 4t 4 [ Cos{2k x )Cos(k y ) 
+ Cos(k x )Cos(2k y ) } + At b Cos(2k x )Cos(2k y ) 



(2) 



and in the chain 



-fj, c — 2hiCos(k y ) + 2h 2 Cos(2k y 



(3) 



The resulting FS, shown in Fig. 1, very closely resembles the observed FS|| and the posi- 
tion of the van Hove singularities (about lOmeV away from the FS|]l7|) seen in Y Ba2Cuj > 0^ 1 
with the following choice of the hopping parameters , 

[li p , t u t 5 ] = [-290, 300, 100, 25, 25, 20}meV and [/x c , h u h 2 ] = [-740, 530, 135}meV. 

The non-interacting Hamiltonian Hq can be easily diagonalised by the introduction of 
new quasiparticle operators «k and (3-^ defined by 



Cka\ v / t 2 + (e+- e )2 v ^ t 2 + (£ -_ £ , ) 2 /a k \ 

V v* 2 +(^ + -) 2 v / * 2 +(— ') 2 / ' 

The quasiparticles corresponding to a kcr (/3 ko .) live on the upper (lower) FS given by e + = 



0(e~ = 0), where e = | (e + e') ± y (e — e') 2 + At 2 denote the dispersions of the two 
quasiparticle bands and the k-dependences have been omitted for notational convenience. 
Fig. 1 shows the corresponding FS in the first quadrant of the Brillouin zone (BZ) and Fig. 2, 
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the partial and total density of states (DOS) with the van Hove singularities predominantly 
located about lOmeV below the FS[|I~7|j. 

Superconductivity occurs through the interaction part of the Hamiltonian fHf 



Hi = ~gJ2 c k't c -k'i c -k'| c kT - g'J2 4'T rf -k'i rf -k'|4t 

kk' kk' 

+K (4 P T rf -k'j. c -k'|C kT + h.c) (5) 

kk' 

where g and g' are the pairing interactions (attractive, non-retarded) in the plane and chain 
respectively; K is the measure of pairing interaction between plane and chain and is assumed 
repulsive. It is this repulsive interaction that forces the quasiparticles in planes and chains to 
stay away in real space forcing J2k( c kad-]<.~a) — which requires the pairing wave function 
to have nearly equal regions of positive and negative signs in the BZ. 
A simple BCS mean-field theory gives the following Hamiltonian, 



Hf< = £ [Act t cL kl + A'4Al + H + + A "^j^ < 6 > 

where A, A' are defined by, 

A = -gJ2(c- H c^} + if^(cL ki d kT ) 

k k 

A' = K^(c_ ki c kT ) + g'J2(d-Hd^) 
k k 

Transforming to the quasi-particle operators, H™ can be re-written as 

H? f = Y, A ^al kl +Y,KPlP j - H + h.c.)+ const. (7) 
k k 

(8) 

Here, A k = A * e +^ e ^7^ anc ^ ^k = A ( £ +!_^)a+7^ • ^ ne ^ wo bands are completely decoupled 
above leading to an identification of A and A' as the gap functions in the plane and the 
chain respectively. On the FS ee' = t 2 and one has a single OP given by A k = A^ = A *+^ e ■ 
Since e, e' have the same sign on the FS, the only way A k changes sign (and consequently 
have a node) along the way on the FS is by having A and A' out of phase. This is, of course, 
ensured by the repulsive interaction K between the two bands as discussed above. 

Diagonalization of the full Hamiltonian Hq + H™-^ is achieved using the usual Nambu 
representation i/j^ = (c^ c_ k | d-ku and writing the mean-field Hamiltonian as ^Adijjk, 
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where M is the 4x4 mean-field Hamiltonian matrix. The four eigenvalues are obtained 
from 



2E± 2 = rj + rj'± y/( v - t/) 2 + 4(r 2 + S 2 ) (9) 

where, rj = e 2 + A 2 + t 2 , rf = e' 2 + A' 2 + t 2 , r = t(e + e') and 5 = t(A - A') 

On minimization of the mean-field free energy with respect to A* and A'* the gap equa- 
tions are obtained as (with (3 as the inverse temperature). 



KA + gA' ^fdE+^__ L {3E+ , dE£ X ___JE V 



and 



99' - K 2 



KA' + g'A 
99' ~ K 2 



? I 5A k ' 



tanh- 



+ 



8A> 



-tank- 



(10) 



22 I -^rr-tanh 



dA 



+ 



<9A 



(11) 



These gap equations are solved numerically and the two gap functions, have the usual 
square root dependence on temperature, but opposite signs. This change of sign of the gap 
function on the same sheet of the FS has important bearing on the physics of the model. 

Since the FS is a combination of contributions from the chain and plane bands, as we 
move along any one of sheets of the FS the OP changes from A to A' or vice versa, and 
hence changes sign. The appearence of OP with different signs can be seen by expanding 
the Free energy, obtained from MF Hamiltonian (Eqn. 6), about (A, A') =(0,0) close to the 
transition. The expression to second order in A and A' is, 



F = a\A\ 2 + b\A'\ 2 + 4c(AA* + A* A) 



(12) 



with a = , g „ 2 — u — v — w,b = , 9 „, — u — v + w;,andc 
where, 



K 



.'/.'/ 



-K 2 



+ V 



4ir 



tanhpE + (0,0) tanh(3E~(0, 0) 



2 k v /(e 2 -e' 2 ) 2 + 4t 2 (e + e')^ 
1 



w 



y 

4 k v /(e 2 -e' 2 ) 2 + 4t 2 (e + e')^ 



£+(0,0) ^"(0,0) 

~tanh(3E+ (0,0) tanhpE~ (0,0) 
£+(0,0) £"(0,0) 

tanhf3E + (0,0) tanh(3E~ (0,0) 



£+(0,0) 



£-(0,0) 
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A nonzero c ensures that there is a single T c , because whenever one of the order param- 
eters approaches zero, it forces the other to vanish simultanously. Since v can be seen to 
be positive and the interaction strength K is not large enough (K 2 < gg'), the sign of the 
coefficient of the third term in the expanded free energy above is always nonzero (unless 
t — K — ) and positive. The free energy is then minimized with the order parameters A 
and A' having opposite signs. This, of course, ensures that there is a node on the FS and 
that there is a single transition temperature. Experimentally, one does indeed see only one 
specific heat anomaly ||. 

Having obtained A and A' thus, it is natural to look for the sign of the OP on the 
FS. We show in Fig. 3 how Ak behaves as one moves along the (non- interacting) FS. The 
consequent node positions on the FS are shown in Fig. 1 for optimal values of g, g' and K 
(that fit the T c , see discussion below). The nodes are very close to the (7r,7r) direction as 
seen experimentally for YBa2Cus07, a feature shared by the d-wave OP. We also calculate 
the position of the nodes as a function of doping (inset of Fig. 1). The node positions turn 
out to be not very sensitive to doping x in YBa 2 Cu 3 6+:2 ;, varying within ±10° about 45° 
i.e., (7r,7r) direction. Although there is no experimental data to compare it with presently, 
this remains a testable prediction of the model at hand. 

It is in the fitness of things that at this point we calculate the dependence of the su- 
perconducting transition temperature T c as a function doping x shown in Fig. 4. We have 
used as fitting parameters g,g' and K, while the band parameters have already been fixed 
to fit the ARPES data as discussed earlier. We chose the values of these parameters such 
that g, g' and K are within reasonable limits giving T c at x = 1 around 100°K. The slightly 
higher value (100°K as compared to the experimental value ~ 90°) that we have chosen is to 
account for the strong fluctuations that are present in these systems owing to their intrinsic 
low dimensionality [|l^]. A reduction of T c by about 10 to 20 percent due to fluctuations is 
generally observed. The optimal values, thus fixed, turned out to be 225, 140 and 125 meV, 
with t = 50meV fixed earlier. The T c is maximum at x — 1 where the Fermi level is closest 
to the van Hove singularity (Fig. 2). As one moves away from this filling, the DOS at the 
Fermi level drops drastically and then moves over a plateau. The T c versus x curve follows 
this pattern as expected. We have also studied the behaviour of T c as a function of the 
interaction parameters g, g' and K and the interband hopping strength t. The variations of 
T c with g and g' show the usual BCS behaviour. Interesting variations are observed in the 
dependence of T c on K and t. As the two gap parameters A and A' have opposite signs, 
and t induces band mixing, T c goes down quite dramatically with t. Conversely, K works 
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counter to this and so pushes the T c up showing up as a minimum (that moves towards lower 
temperature as t is brought down from 50 to 15meV) followed by the usual rise as expected 
in a two band model for superconductivity |T^]. 

The two different order parameters in our model are not related by any symmetry, except 
that they have different signs. The c-axis Josephson current, therefore, does not cancel out 
completely as in the case of a d-wave order parameter, although partial cancellation occurs 
due to the opposite signs of A and A' resulting in a small (compared to an s-wave OP) 
but finite current as seen in experiments ||. We found the node on the FS close to the 
(•7T, vr) direction, which is not very sensitive to the doping concentration. The appearance of 
the node on the FS ensures a linear temperature dependence of penetration depth and low 
energy excitations seen in the thermodynamic measurements. The order parameter in a— ( 
b— ) direction is mainly chain (plane) like and hence provide a natural explanation for the 
observed phase change rr in SQUID experiments. 

The s-wave order parameter is weakly influenced by impurities (Anderson theorem). 
In a d— wave superconductor, however, scattering between lobes of different sign on the 
FS reduces the T c much faster. This averaging out of the sign and consequent sensitivity 
towards non-magnetic impurities will occur in the present model as well, albeit at a much 
slower rate than in the d— wave case. Sign of OP changes in the present model only by 
scattering from planes to chains, involving strong scattering processes only. Sensitivity of 
T c towards non-magnetic impurity will, therefore, be intermediate between an s— and a d— 



wave superconductor in our case ||20|| , as seen in the high T c materials ^, IT]. The superfluid 
density, obtained by calculating the leading order fluctuations over the mean-field state gives 
an anisotropic penetration depth ||20|| . 

In the YBCO system, the unit cell contains two planes and one chain. So that a more 
realistic model would be where these additional features are present. It is important to note, 
however, that the conclusions of the present model are not very sensitive to whether we work 
with a chain and a plane, two or more planes or chains, or a combination of them. This a 
simple model that incorporates the effect of having more than one OP residing on different 
bands that are coupled and shows that with the correct choice of phase of the different OP 
can capture many of the unusual features of the symmetry of OP in these systems. It would 
be nice to see experimental results on the position of the node(s) on the FS as a function 
of doping. In the d— wave case, since the position is dictated by symmetry, there will be no 
change as the FS changes shape, while in any of the s— wave models, it will. Though in our 
model, for a range of doping, the sensitivity is not too strong. 
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In conclusion, we have shown that it is possible to understand many of the contradictory 
experimental evidences from a simple and realistic model applicable to layered superconduc- 
tors. 
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FIGURE CAPTIONS 



Fig. 1 The two sheets of FS are shown in the first quadrant of the BZThe dotted(solid) 
line corresponds to e~ = (e + = 0). The filled circles show the position of nodes of OP on 
FSat x = 0.9. The variation of the node positions on FS with x is shown in the inset. 

Fig. 2 The DOS for noninteracting quasiparticles is shown. vHS is around 10 meV below the 
Fermi Energy. The inset shows the partialal DOS, where the dotted(solid) line corresponds 
to the (T ( e + ) band. 

Fig. 3 The anisotropic OP comes out to be same in two different sheets of FS.In the figure 
the dashed (solid) line shows how A fc (= A k ' behaves on FS corresponds to e~ = 0(e + = ) 
at x = 0.9. Clearly two point nodes are visible. 

Fig. 4 The behaviour of T c with doping concentration, shown in this figure is somewhat 
similar to the experimental curve. The platue around x = 0.8 is prominant and the max. T c 
about x 1 is because of the presence of vHS there. 
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